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Abstract 



The interlayer pair tunneling model of Anderson et al. is generalized to 
include the strong coupling effects associated with in-plane interactions. The 
equations for the superconducting transition temperature T c are solved nu- 
merically for several models of electron-optical phonon coupling. The non- 
magnetic in-plane impurity scattering suppresses T c in all cases considered, 
and it is possible to obtain a fair agreement with experiments for a reasonable 
choice of parameters. For the anisotropic electron-phonon coupling proposed 
by Song and Annett we find that the interlayer pair tunneling can stabilize 
the d;j.2_ y 2-wave superconducting state with a high T c . Moreover, in this case 
there is a possibility of an impurity induced crossover from the d a ,2_j / 2-wave 
state stabilized by the interlayer tunneling to the s-wave state at a low impu- 
rity concentration. We also calculate the isotope effect associated with the in- 
plane oxygen optic mode and its dependence on the strength of the interlayer 
pair tunneling. Small positive values of the isotope exponent are obtained for 

strengths of pair tunneling that give high transition temperatures. 
PACS numbers: 74.20.-z, 74.20.Mn, 74.72.-h 
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I. INTRODUCTION 



One of the theories for high-T c copper-oxide superconductors is the interlayer pair tun- 
neling (ILPT) model first proposed by Wheatly, Hsu and Anderson |l|] and later refined by 
Chakravarty, Sudb0, Anderson and Strong [0-0. Recently, Chakravarty and Anderson |3j 
gave an indirect justification of this model based on a non-Fermi liquid form of the normal 
state electron propagator. In the ILPT-model the pairing in the individual copper oxygen 
layers is enhanced and sustained by the pair tunneling between the layers within the unit 
cell. The specific symmetry of the component of the order parameter resulting from in-plane 
interactions is not an essential feature of the model, although in the original work || it was 
assumed that this component has s-wave symmetry. In this paper we generalize the BCS- 
like form of the ILPT-model given by Chakravarty et al. to include the retardation (i. e. the 
strong coupling) effects resulting from in-plane interactions. This generalization is neces- 
sary in order to obtain a more realistic dependence of T c and other quantities characterizing 
the superconducting state on the interaction parameters 0. Moreover, the strong coupling 
form of the interlayer pair tunneling model is suitable for including the effect of in-plane 
nonmagnetic impurity scattering. The dependence of T c on impurity concentration 
is considered to be an important indicator of the symmetry of the order parameter in oxide 
superconductors and of the underlying pairing mechanism [0,0 • 

The self-energy equations derived in this paper are valid for any kind of in-plane pairing 
interaction within the one-boson exchange approximation. In our numerical work, how- 
ever, we consider only the case when the in-plane pairing arises from electron coupling 
to optical phonons. This was motivated by the fact that the ILPT-mechanism is novel 
enough that its consequences should be examined first when the in-plane pairing is caused by 
the conventional electron-phonon interaction before more exotic in-plane interaction mod- 
els are considered. Also, Song and Annett |14]] recently derived an effective single-band 
Hubbard-type Hamiltonian for Cu02 planes which includes the electron-phonon coupling to 
oxygen breathing modes. The electron-phonon matrix element squared is proportional to 
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sm 2 ((k x — k' x )/2) + sm 2 ((k y — k' y )/2), where k and k' are the electron momenta. With this 
form of coupling Song and Annett initially predicted that the order parameter with d x 2_ y 2- 
wave symmetry would lead to a higher transition temperature than the order parameter 
with s-wave symmetry, because in the former case the on-site Coulomb repulsion U becomes 
ineffective. Subsequently, they found this prediction to be erroneous, which we indepen- 
dently confirmed during the course of this study. We found, however, that the d x 2_ y 2-wave 
state could be stabilized by the interlayer tunneling. 

The rest of the paper is organized as follows. In Section II we give in some detail 
the derivation of the interlayer tunneling contribution to the electron self-energy in the 
superconducting state and list the well known results for the self-energy parts arising from 
in-plane interactions. In addition, we summarize the T c -equations for the case of pairing 
induced by electron-optic phonon coupling. Section III contains our numerical results for the 
transition temperature as a function of increasing disorder and the oxygen isotope exponent, 
and finally in Section IV we give conclusions. 



II. THEORY 

A. Self-energy due to interlayer pair tunneling 

In the ILPT-model it is assumed that in the superconducting state the quasiparticle 
picture is approximately valid for motion within a layer, while the coherent motion of quasi- 
particles from layer to layer within the unit cell is blocked [@,|J. The part of the Hamiltonian 
which describes the interlayer pair tunneling for two layers per unit cell [^],[5| is 

Hj = ~E TjM [c^c%c%^ + H.c] , (1) 

k 

where offi is an electron creation operator for the state of two-dimensional momentum k 
and spin j in the layer (z), and 

t 2 

Tj(k) = — [cosf^a) - cos(k y a)] 2 , (2) 



as suggested by band structure calculations JT^|. Here, t± characterizes the high energy 
single-electron coherent hopping from layer to layer, and is estimated to be between 0.1 eV 
and 0.15 eV ||. The parameter t enters the tight-binding dispersion for the electron motion 
within a layer 



= —2t[cos(k x a) + cos^a)] — 4t'cos(A; 3 



a) cos^a) — fi 



(3) 



where /i is the chemical potential and a is the lattice constant. To find the contribution to 
the anomalous electron self-energy from the Hamiltonian (1), it is convenient to rewrite Hj 



in the Nambu formalism fl6fl . Introducing the Nambu fields 
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the Hamiltonian (1) could be written as 
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where fi and r2 are the two off-diagonal Pauli matrices ||16|| . This expression looks like the 
sum of two two-body interaction Hamiltonians with the interaction line — Tj (k)5k,k'^q,o an d 
with the interaction vertices fi and T2, respectively []16| |. Since the ILPT-model does not 
consider correlations of the type — (T r (\E'^(r) 1 I / ^ 2 ^(0))), where T T is the Wick's time-ordering 
operator, we consider only the contribution to the electron Nambu self-energy arising from 
the Hartree-type diagrams shown in Fig. 1. The contribution of these two diagrams to the 
irreducible Nambu electron self-energy in layer (1) is 



i 1} (k) 



(5) 



where T is the temperature in energy units, G^(k, iui m ) is the Nambu 2x2 electron Green's 
function for layer (2) at wave vector k and the fermion Matsubara frequency iw m , Tr{- ■ •} 
is the trace, and the overall minus sign arises from one closed fermion loop. With the usual 
form for the total irreducible electron Nambu self-energy in layer i 
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(k, iu n ) = ico n (l - Z« (k, iw n )) f + 0» (k, ^ n )f! + 0« (k, io; n )f 2 + X (i) (k, ^ n )f 3 , (6) 

where Z^ 1 ' is the renormalization function, <jy^ and <^ are the real and imaginary parts, 
respectively, of the pairing self-energy, and is the part of diagonal self-energy which is 
even in iu n , Eq. (5) takes the form 

E«(k) = Tj(k)Tj2 ^(k,^)f 1+ ^)(k,^ 

™ {u m Z®Qs.,iu m )y + (4 2) + X (2) (k,^ m )) 2 + |0( 2 )(k,^ m ) +#(2)(k,^ m )P 

(7) 

It should be noted that the interlayer pair tunneling does not lead to any frequency de- 
pendence in the self-energy, but contributes directly only to the pairing (i. e. off-diagonal) 
self-energy and, as emphasized in || , the resulting self-energy is local in k. Moreover, in the 
weak coupling approximation for in-plane interactions = 1, x — and and 
do not depend on Matsubara frequency, and the sum in (7) could be easily performed using 
contour integration | lq| . One finds 



E«(k) = Tj{k) (^(k)f! + (2) (k)f 2 ) ^tanh (^) , (8) 



where Ev = 

V( e k } ) + |# 2) (k) +^ (2) (k)| 2 , which is the same result (with </>( 2 ) gauged 
away) as that obtained by Chakravarty, Sudb0, Anderson and Strong [Q. 

B. Self-energy due to in-plane interactions 

The precise form of the self-energy due to in-plane interactions depends on the model 
used and we will restrict ourselves to the case where the pairing interaction is due to one- 
boson (e. g. phonon or spin-fluctuation) exchange. The electron-phonon contribution to the 
self-energy is |16|,[7|,|T7 



T 

Eg(k,io; n ) = - — Yl l#k-k',A| 2 £>A(k - k',iu n - iu m )f 3 G {l) (k' , ioo m )f 3 , (9) 

k',A,m 



where N 2 is the number of lattice sites, (?k-k',A is the electron-phonon matrix element for 
the momentum transfer k — k' and the phonon polarization A, D\(k — k', iv m ) is the corre- 
sponding phonon propagator at the boson Matsubara frequency iv m , and is the electron 



propagator for the layer i. An analogous expression is obtained for the self-energy due to 
the exchange of antiferromagnetic spin-fluctuations, except that \g\ 2 D is replaced by the 
spin-fluctuation propagator and the Pauli matrix T3 by the unit 2x2 matrix . 

In the case of phonon-mediated superconductivity one can include the effect of short- 
range Coulomb repulsion within an effective single band Hubbard model for copper-oxygen 



planes ||14j| . The resulting contribution to the electron self-energy is 



T 



^ ) = -^E^GS(k , ,^ m )f 3 



N 2 



(10) 



k',m 



where U is the on-site Coulomb repulsion and is the off-diagonal part of the electron 



Nambu Green's function 17 



Finally we consider the effect of in-plane electron- impurity scattering. We will not con- 
sider all the possible effects of electron-impurity scattering in two-dimensional superconduc- 
tors (e. g. the enhancement of the Coulomb repulsion [B5|), but will confine ourselves to the 
simplest treatment using either the second Born approximation or the t-matrix approxima- 
tion. In the second Born approximation and assuming a constant electron-impurity matrix 
element Vi = (k|V/v|k'), where Vn is the change in the crystal potential due to nonmagnetic 
impurity, the electron self-energy resulting from scattering off impurities is 

n f y? 



N 2 



(11) 



while in the t-matrix approximation the self-energy is given by 



(12) 



C. T c -equations for the case of pairing due to electron-optic phonon coupling 

We consider the case where the in-plane pairing is mediated by an optic phonon of energy 
Qe- Near T c the pairing self-energies become infinitesimal and the self-energy equations could 
be linearized as described, for example, in [|l7],|n|. Assuming that the electron self-energies 
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in each of the two layers are identical (e. g. 0^(k, iu n ) = 0^(k, iu n ) = 0(k, iu n )) and 
defining 

u(k,n)= . (13) 
yj {uj n Z(k, iu n )) 2 + (e k + x(k, iu n )) 2 

the T c -equation reduces to the eigenvalue problem of a real symmetric matrix 

u(k, n) = K(K n; k', m)u(k', m) . (14) 

k',m 

The matrix K consists of several parts associated with various interactions 

K = K ep + Kj + K c + K t . (15) 
The electron-phonon contribution is 

K ep (k,n;k',m) = -j^J2\gk,k',x\ 2 ~ ( _ ^f 2 T n2 £(k, n)ff(k', m) , (16) 

where 

5(k, n) = 1 . (17) 

y (w„Z(k, «<;„)) + (e k + x(k, «<;„)) 

The interlayer pair tunneling contribution is 

Kj(k, n; k', m) = 5^TTj(k')S(k, n)S(k', m) , (18) 

while the contribution due to on-site Coulomb repulsion is 

K c (k, n; k', m) = ~US(^, n)S(k', m) . (19) 

Treating the in-plane impurity scattering in the second Born approximation gives 

n-V 2 

Ki(k, n; k', m) = ^5 n , m S(k, n)S(k', m) , (20) 

while the t-matrix approximation gives 

n-V 2 

K t (k, n; k', m) = -j^5 n , m S(k, n)S(k', m)/D(n) , (21) 



where 
Din) 
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At a given temperature the functions Z(k,iu n ), x(k, io; n ) and the chemical potential /x (see 
Eq. (3)) are determined self-consistently by solving a set of equations 



Z(k, iu n ) = 1 + Z ep (k, icu n ) + Zi(k, ioj n ) 
X(k, za; n ) = x ep (k, zu; n ) + Xi(k, iw n ) 



(23) 
(24) 



together with the equation representing the particle number conservation [20 
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(25) 



where n is the band filling factor and is the (1, l)-component of the electron Nambu 
Green's function. In Eqs. (23-24) Z ep and Xe P are given by 

^k< A 



Z ep (k, iu n ) 



1 UJ n — (jj m 



i 2 + ni 



■u m Z(k, iuj m )S(k',iuj m y 



Xe P (k, iu r , 



2Q, 



T 



e k ' + x(k',^ m )) ^(k^m) 2 



(26) 
(27) 



while Zj and Xi are given by 



Zj(k,za; n ) 



i ' i 



J2Z(k',iuj n )S(k',nf 
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in the second Born approximation, and by 
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AT2 



ft 



(31) 



q \ " q 

in the t-matrix approximation. The transition temperature T c is determined as the highest 



temperature at which the largest eigenvalue of matrix K is equal to 1 (see Eq. (14)). 



III. NUMERICAL RESULTS 



In the numerical calculations we have taken (for definiteness) the same band parameters 
as those in the work of Chakravarty et aL; namely t = 0.25 eV and t'/t = —0.45. The band 
filling factor was set at n = 0.375 corresponding to 0.75 electrons per cell. In Fig. 2 we 
show the density of states N(E) for these parameters obtained by adapting the tetrahedron 
method to a 400 x 400 square lattice. We assumed that the electrons couple to an optic 



phonon of energy = 62 meV (i. e. 500 cm ) 13]. Two models for the electron-phonon 



matrix element were considered. In the first model, which we will refer to as the isotropic 



2 

model, the momentum-independent |<7k,k/,A| was assumed 

X!l£k,k',A| 2 = M 2 , (32) 

A 

2 

and the value of \g\ was chosen such that the electron-phonon mass renormalization pa- 
rameter A ~ 0.5 . In the second model, which we will refer to as the anisotropic model, the 
momentum dependence of |gk,k',A| 2 was taken to have the form given by Song and Annett 
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(33) 



with the same value of \g\ as in (32), so that the maximum in (33) is equal to \g\ in the 
isotropic model. Because the interlayer pair tunneling contribution to the kernel in the In- 
equation (14) is local in k, Eq. (19), the calculations had to be performed in k-space except 
when considering isotropic in-plane interaction with Tj = 0. In this case it is possible to 
convert the k-sums into integrals over the electron energies and use the electronic density of 
states calculated for a large (400 x 400) lattice. The results for Tj = and isotropic in-plane 
interaction served as a check of the accuracy of the results obtained from the calculations 
in k-space. Due to memory size restrictions on the computers that were available to us (4 
processor SGI R4400 and Fujitsu VPX240/10) the largest lattice size that we could consider 
was 64 x 64. We found that it is absolutely critical to add and subtract the noninteracting 
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form of the band filling factor n to the expression Eq. (25) and to evaluate the added part as 
/ dEN(E) I (exp((E — /i)/T) + 1). Otherwise, the truncation of the sum over the Matsubara 
frequencies in (25) and the finite lattice size could lead to an error in T c as great as 55% in 
the case of isotropic in-plane interaction with Tj = 0. This error in T c is largely due to error 
in the chemical potential fj, which leads to an incorrect value for the density of states near 
//. We found that this trick of adding and subtracting leads to values of T c that are accurate 
to better than 5% for the largest lattice size that we could consider. The largest eigenvalue 
of the matrix K in (14) was obtained using the power method [^] and due to the simple 
structure of the sums over the Matsubara frequencies in (25-31) there was no need to use 
the fast Fourier transform technique of Serene and Hess p3| . The resulting code vectorized 
93-97% on Fujitsu VPX240/10. 

A. T c suppression by in-plane impurity scattering 

We first consider the isotropic model of electron-phonon interaction and, naturally, as- 
sume the s-wave symmetry of the pairing self-energy. Fig. 3 illustrates the suppression of T c 
by in-plane impurity scattering obtained within the Born approximation for four different 
values of the interlayer pair tunneling parameter Tj = t\jt (Eq. (2)), and with the on-site 
Coulomb repulsion U = 0. It should be stated from the outset that with the strong cou- 
pling effects (i. e. renormalization) one needs a larger value of Tj to achieve a transition 
temperature of about 100i^ with no disorder than in the BCS-like treatment PJ24]] (here the 
electron phonon mass renormalization parameter is A = 0.48 as deduced from the value of 
Z at the first Matsubara frequency and Qe — 62). In the Born approximation the impurity 
scattering is parametrized by niVj 2 , and we plot T c as a function of this quantity. If we take 
l/2r; = 7cN(Ep)niV^, where N(Ep) is the band electronic density of states at the Fermi 
level, as the measure of the elastic scattering rate, the range shown in Fig. 3 corresponds to 
about HOmeV"; with Vi . = t = 250meV the maximum value of niV^ in Fig. 3 is obtained for 
the in-plane impurity concentration = 0.48 per cell. The overall shapes of the curves in 
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Fig. 3 are similar to the results obtained by Bang |24| in the BCS-type of treatment using 
the circular Fermi surface and 7j(k) oc | cos20|, where <fi gives the position of k on the 
Fermi surface. However, we find that T c is suppressed at a much slower rate than what was 



obtained by Bang in the Born limit [p^j . At first, it is surprising that we get a drop in T c 
with increasing riiV- 2 for Tj = 0. In this case there is no gap anisotropy which could be 
washed out by the impurity scattering leading to the suppression of T c . Also, the structure 
in N(E) within a range ±Qe around the Fermi level does not seem to be significant enough 
(see the inset in Fig. 2) for the smearing caused by the elastic scattering rate l/2rj = 0(Qe) 
to have any significant effect on T c . We checked the result for Tj = by converting the 
k-sums into integrals over electronic energies, as discussed at the beginning of this section, 
and found the same result. Upon inspection we found that with increasing n^Vf there is 
a slight shift in the chemical potential to the region of lower density of states. Although 
the reduction of the density of states at the chemical potential is small, the exponential 
dependence of T c on the interaction parameters presumably leads to the observed decrease 
in T c for Tj = 0. Note that the rate of suppression is greater for increased Tj. The reason 



for this trend has been discussed by Bang 



Next we consider the effect of in-plane disorder in the t-matrix approximation. The 
results are shown in Fig. 4 for Tj = 90meV. We plot T c as a function of the in-plane 
impurity concentration rii for several values of the impurity scattering potential parameter 
Vi. Note that for V{ = t = 250meV the t-matrix approximation (solid line) and the Born 
approximation (dots) give very similar results as one would expect in the limit of small Vi 
(see Eqs. (11-12)). Increasing Vi leads to a more rapid suppression of T c with increasing 
impurity concentration, and the unitary limit is reached by V, = 6t-8t with U = 0. Note, 
however, the change in curvature of T c versus rij as Vi is increased. This trend was not 
found in the weak coupling calculation of Bang J^] in crossover from the Born limit to the 



unitary limit. However, the overall rate of suppression of the transition temperature with 
increasing disorder that we find in the unitary limit is comparable to the rate found by 



Bang |24] for A = 0.5 (with our N(Ep) = 1.16 x 10 3 states/cell/meV/spin the impurity 
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concentration ni = 0.1 cell -1 corresponds to the impurity scattering rate in the unitary limit 
Tj = rii/ir /N(Ep) = 27 meV). We also obtained step-like features in the T c -curves in the 
unitary limit which we are not able to associate with any particular feature of the model 
and/or the numerical procedure used. The experiments p|-|TTf in general do not report the 
data on the very fine scale over which we observe the steps, and only in |J was there an 
attempt to interpret fine features of the observed dependence of T c in Yi^Pr^E^Cu/iOs on 
Pr concentration. It should be kept in mind that only the experiment of Tolpygo et al. |TT|] 



addresses specifically the in-plane defects in YBa2Cu306+z at the fixed carrier concentration 
to which our model calculations apply. The most important aspect of Fig. 4 is that it 
illustrates the profound effect of the Coulomb interaction on the dependence of T c on the 
concentration of in-plane impurities. For U = 8t (the bandwidth) the solid curve in Fig. 
4 obtained for Vi = t in the t-matrix approximation is pushed down to the line given by 
squares. The decrease in T c for U = 8t is about 5K per 1% of in-plane defects, similar to 
the value found by Monthoux and Pines ]nj for Vi = t in the model of spin-fluctuation- 
induced superconductivity and d-wave paring, and to the value measured by Tolpygo et al. 
\T1\. Moreover, we found that for this choice of parameters (A = 0.48, U = 8t = 2eV) 



switching off the interlayer pair tunneling reduces the transition temperature from 84. 5K 
(for Tj = 90meV) to 1.6K for rij = 0. This illustrates the remarkable effect of the interlayer 
pair tunneling mechanism on the enhancement of T c . 

Next, we turn to the anisotropic model of the electron-phonon coupling function, Eq. 
(33). The results for Tj = 90meV using the t-matrix approximation are shown in Fig. 5. 
As we have mentioned in the Introduction we were not able to obtain a finite transition 
temperature assuming d x 2_ y 2 -symmetry of the pairing self-energy for Tj = down to the 
lowest temperature that we could consider using the k-space method (about 20K). However, 
with Tj = 90meV we obtained a transition temperature of 114K assuming d x 2 .^-symmetry 
of the gap for rti = 0. It is interesting that the s-wave case with U = At and the d x 2_ 3/ 2-wave 
case (the on-site Coulomb repulsion drops out, Eq. (10)) give quite a similar dependence of 
T c on rii for both Vi — t and Vi = lOt. The results obtained for Vi = t are similar to the 
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experimental results on YBa2Cus06+z with in-plane oxygen defects ||11|| , although we find 
that the squares in Fig. 4 more closely resemble the experimental data at the highest values 
of Hi where the data seem to fall on a curve that becomes less steep as rii is increased. In the 
unitary limit Vi = lOt the T c -curves initially rise with increasing rti and then precipitously 
drop. There seems to be a common threshold 7ij beyond which superconductivity disappears 
for both the s-wave pairing with either U = At or U = 8t and for d x 2_^2-wave pairing. We 
have found a similar behavior for s-wave pairing with U = 8t, Tj = 90meV and V, = t (not 
shown here), except that the initial rise in T c is much less pronounced and the threshold 
occurs at a higher value of rij. 

We would like to point out that with the electron-phonon coupling function given by 
Eq. (33) it is possible to have an impurity induced crossover from the d^^-wave state in 
a very pure system to the s-wave state at a higher impurity concentration, Fig. 6. All of 
the interaction parameters for the two curves in Fig. 6 are the same. The only difference is 
that for the solid curve the d x 2_ y 2-symmetry of the pairing self-energy is assumed (in this 
case the on-site Coulomb interaction drops out, Eq. (10)), while for the triangles the s-wave 
symmetry is assumed. At a given impurity concentration the system will go into a state 
with a higher T c in order to lower its free-energy. 



B. The isotope effect associated with the in-plane oxygen optic mode 

We have examined the isotope effect associated with the optic phonon which mediates 
the in-plane interaction. In the original work of Chakravarty et al. it was suggested that the 
interlayer pair tunneling mechanism could explain a small isotope effect in high-T c copper 
oxide superconductors simply because in the interlayer tunneling model the most important 
pairing process is associated with the pair tunneling. Our results for the isotope exponent 
a = —d\nT c /dhiM associated with the optic mode at VLe are shown in Fig. 7. In the same 
figure we give the corresponding transition temperatures. The impurity scattering was set 
equal to zero and the results for s-wave symmetry were obtained for the on-site Coulomb 
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repulsion U equal to zero. Note that for the isotropic model of electron-phonon interaction, 
Eq. (32), we get the classical result a = 0.5 for Tj = 0. In the same model Tj = 90meV 
gives T c = 108K and a = 0.18. Turning on the on-site Coulomb repulsion to U = 8t (the 
bandwidth) reduces the transition temperature to T c = 84. 5K and the isotope exponent to 
a = 0.05 (not shown in Fig. 7)-a value approximately equal to what is found for the oxygen 
isotope effect in high-T c Y-Ba-Cu-0 systems p5| |. It should be mentioned, however, that 



in the site selective oxygen isotope experiments of Nickel et al. [Ej| , where only the oxygen 
in copper-oxygen planes is replaced by heavier isotope, a small negative isotope effect is 
observed with the partial isotope exponent a = —0.01 ± 0.004-close to the resolution limit 



[ 25 1 . For anisotropic electron-phonon interaction, Eq. (33), and assuming s-wave symmetry 
of the gap we generally get lower values of a than in the isotropic case. For Tj = 90meV and 
U = the transition temperature is 123K and a = 0.04. Turning on U = 8t reduces the T c 
to 78K and increases the isotope exponent to a = 0.08. For d^.^-symmetry of the gap we 



obtain very small positive values of a-likely smaller than the experimental resolution [EH] . 



IV. CONCLUSIONS 

We have generalized the interlayer pair tunneling model of Anderson and coworkers 
to include the retardation effects associated with in-plane interactions. Through numerical 
solutions of the T c -equations for a model in which electrons couple to an optic phonon at 500 
cm -1 (i. e. 62 meV) we found, without trying to fit the experiments, that a reasonable choice 
for the band parameters (t = 250meV, t'/t = —0.45), band filling factor (0.75 electrons per 
cell), electron-phonon coupling (A = 0.48), on-site Coulomb repulsion (U ~ the bandwidth), 
and the interlayer pair tunneling strength (t± = 0.15eV) yields results in surprisingly good 
agreement with the experiments on both a T c suppression by in-plane oxygen defects [T1J and 
the oxygen isotope effect [E_5j in YBa2Cus06+x- The best agreement is found for the isotropic 
model of the electron-phonon coupling function with U = 8t which leads to a T c suppression 
rate of about 5K per 1% of the in-plane defects (with the impurity matrix element Vi = t) 
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and to the oxygen isotope exponent a = 0.05. This case also best illustrates the importance 
of the interlayer pair tunneling process in raising the transition temperature, since reducing 
t± from 0.15eV (i. e. Tj = 90meV) to zero decreases the T c from 84. 5K to 1.6K. We also 
found that for the anisotropic form of the electron-phonon coupling proposed by Song and 
Annett Wffl, Eq. (33), the interlayer pair tunneling can stabilize the superconducting state 



with d x 2_ ?/ 2-symmetry at a high T c . This stabilization occurs because the pair tunneling 
contribution to the pairing self-energy is local in k. Morever, it is possible to have impurity 
induced crossover from the d x 2_ y 2-state in a "perfect" sample to the s-wave state at a higher 
concentration of defects. This is illustrated in Fig. 6 for Tj = 90meV and U equal to half 
the bandwidth, but we have also found examples of such a crossover for other values of Tj 
and U. 
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FIGURES 
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FIG. 1. The inter layer pair tunneling contribution to the electron Nambu self-energy. 
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FIG. 2. The electronic band density of states for the dispersion given by Eq. (3) calculated 
for 400 x 400 lattice using the tetrahedron method. The energy is measured in units of t and 
Ep indicates the Fermi level. The inset shows the density of states (using the same units) in the 
interval [Ep — ^e, Ep + ^e], where Qe = 62meV is the energy of the optic phonon. 
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FIG. 3. The dependence of T c on n^ 2 , where i%i is the concentration of in-plane impurities 
and Vi is the electron-impurity matrix element, calculated in the Born approximation, Eq. (11), 
for several values of the interlayer pair tunneling strength Tj = t\/t, assuming the isotropic 
electron-phonon coupling model (32) (A = 0.48), s-wave pairing, and the on-site Coulomb repulsion 
U = 0. 
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FIG. 4. The dependence of T c on the concentration rij of in-plane impurities for the isotropic 
electron-phonon coupling model (32) (A = 0.48) and Tj = 90meV calculated in the t-matrix 
approximation, Eq. (12), for several values of the electron-impurity matrix element V*, except for 
the results given by filled circles which were obtained in the Born approximation (the same as the 
solid curve in Fig. 3) . S-wave symmetry is assumed and the on-site Coulomb repulsion U = except 
for the results given by squares where U is set equal to the bandwidth. 
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V i= t, U=4t (s-wave) 

V i= 10t, U=4t (s-wave) 
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FIG. 5. The dependence of T c on the concentration raj of in-plane impurities for Tj = 90meV 
and for the anisotropic electron-phonon coupling model (33) (the same value of \g\ 2 was used as 
in Figs. (3) and (4)) calculated in the t-matrix approximation, Eq. (12), for the electron-impurity 
matrix element V{ set equal to either t or Wt (the unitary limit). The assumed symmetry of the 
gap is indicated in the brackets and the value of the on-site Coulomb repulsion U is indicated in 
the legend. For d x 2_ y 2 -symmetry of the gap U drops out, Eq. (10). 
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FIG. 6. The dependence of T c on the concentration rij of in-plane impurities for Tj = 90meV 
and for the anisotropic electron-phonon coupling model (33) (the same value of \g\ 2 was used as in 
Figs. (3), (4) and (5)) calculated in the t-matrix approximation, Eq. (12), for the electron-impurity 
matrix element V* = t. The results for d x 2_ y 2 -symmetry of the gap are given by the solid line and 
the results for s-wave symmetry are given by triangles. Note that as rii increases the state with 
the highest T c crosses over from d x 2_ y 2 -symmetry to s-wave symmetry. 
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FIG. 7. The isotope coefficient a = —d\a.T c /dhiM associated with the oxygen optic mode at 
500 cm -1 and the corresponding T c as functions of the interlayer pair tunneling strength Tj = t\/t. 
In all cases the on-site Coulomb repulsion U is equal to zero. The results obtained with the isotropic 
electron-phonon coupling model (32) are labeled as (s-wave isotropic). The results obtained with 
the anisotropic electron-phonon coupling model (33) are labeled as (s-wave anisotropic) or as 
(d-wave anisotropic) depending on whether s-wave or d^^-wave symmetry is assumed, respec- 
tively. The value of \g\ 2 is the same as in Figs. (3-6) leading to A = 0.48 in the isotropic case. 
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